clear all;
clc;

Re = 6378.145e+3;

t0 = 0;
tf = 24*3600;
tspan = [t0, tf];

r0 = [0;0;0];
v0 = [0;0;10];
x0 = [r0; v0];

option = odeset('RelTol', 1e-5);
[t_orbit, x_orbit] = ode45('dif_unBallon', tspan, x0, option);
%%%%%%%%%%%%%%%%%%%%%%%%%%
figure
plot3(x_orbit(:,1),x_orbit(:,2), x_orbit(:,3));
figure
plot3(-x_orbit(:,2),x_orbit(:,1),x_orbit(:,3));
figure
plot(-x_orbit(:,2),x_orbit(:,1));
